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X-ray Emission Properties of the Old Pulsar PSR B2224+65 
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Abstract. Using archival Chandra data we studied the X-ray emission properties of PSR B2224+65 and its environment. 
Albeit limited by photon statistics the spectral analysis suggests that the bulk of the emission from PSR B2224+65 is non- 
thermal. Fitting a power-law model to the observed energy spectrum yields a photon index of T = 1 .58^2 m . The possible origin 
of the non-thermal pulsar emission is discussed in the context of the outer-gap model. We did not find any evidence for a 
compact nebula around PSR B2224+65 though the Chandra data reveal the existence of an extended feature which appears to 
be associated with PSR B2224+65. It extends from the pulsar position about 2 arcmin to the north-west. Its orientation deviates 
by ~ 118° from the pulsar's proper motion direction. Investigating its energy spectrum shows that the emission of this extended 
feature is much harder than that of the pulsar itself and is non-thermal in nature. 
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1. Introduction 

High energy radiation properties of rotation-powered pulsars 
are known to approximately correlate with their spin-down age. 
X-rays from young pulsars with age less than ~ 10 5 years are 
believed to be dominated by magnetospheric emission, charac- 
terized by a power-law spectrum and narrow pulse profiles. For 
middle-aged pulsars in the age bracket between ~ 10 5 — 10 6 
years the pulsars' X-radiation is dominated by cooling emis- 
sion from the hot stellar surface. This class of pulsars is ob- 
served to have composite spectra consisting of a soft ther- 
mal component, a harder thermal component from heated po- 
lar caps as well as some non-thermal emission. The thermal 
and non-thermal contributions in the emission of younger and 
middle-aged pulsars has already been well disentangled and 
studied by previous X-ray observatories (for recent reviews see 
e.g. Becker & Aschenbach 2002; Kaspi et al. 2004 and refer- 
ences therein). 

Old neutron stars with ages > 10 6 years have already 
cooled down to less than five hundred thousand degrees so 
that the thermal emission from their surface has faded from 
view for X-ray observatories. Recently, the much improved 
sensitivity of the XMM-Newton and Chandra observatories en- 
abled a series of comprehensive studies to probe and iden- 
tify the origin of the X-radiation from old rotation-powered 
pulsars. Until now, results have been presented for the pul- 
sars B0950+08, B0823+26, J2043+2740 (Becker et al. 2004), 
B0628-28 (Becker et al. 2005), B0943+10 (Zhang, Sanwal & 
Pavlov 2005), Bl 133+16 (Kargaltsev, Pavlov & Garmire 2006) 
and for PSR B1929+10 (Becker et al. 2006). For all old pulsars 
which were observed with sufficient photon statistics the X-ray 
emission turned out to be dominated by non-thermal emission 



originating in the pulsar's magnetosphere. X-ray pulses, if de- 
tected, often show narrow features. 



A very interesting member in this class of old pulsars 
is PSR B2224+65. The proper motion of PSR B2224+65, 
jUra=144 mas/yr and /z Dec =112 mas/yr, is among the highest 
observed by now (Manchester et al. 2005). The pulsar's inter- 
action with the ISM produces a magnificent bow shock nebula 
which was discovered in Ha by Cordes et al. (1993). Because 
the shape of this nebula resembles that of a guitar the nebula 
got dubbed Guitar Nebula. The pulsar has a period of about 
0.68 s and a period derivative of 9.55 x 10~ 15 s/s, implying a 
spin-down luminosity of £ = 1 . 19 x 10 33 erg/s and a spin-down 
age of somewhat more than one million years. The distance to 
PSR B2224+65 is not very well constrained. The radio dis- 
persion measure based distance in the Cordes & Lazio (2002) 
Galactic free electron density model is 1.86 kpc. However, 
Chatterjee & Cordes (2004) found by modeling the head of the 
Guitar Nebula that the pulsar distance could be as close as 1 
kpc. For this work we adopt the distance of 1 kpc as suggested 
by Chatterjee & Cordes (2004). The ephemeris of this pulsar 
as adopted from Manchester et al. (2005) are summarized in 
Table 1. 



In this paper we present a detailed X-ray study of PSR 
B2224+65 and its environment based on archival Chandra data. 
Brief results from this observation were already reported by 
Wong et al. (2003) and Zavlin & Pavlov (2004). The structure 
of this paper is as follows. In §2 we describe the observations 
and data analysis and in §3 we summarize the results which for 
the pulsar are discussed in the context of the outer-gap model. 
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Fig. 1. Chandra ACIS-S3 image in the energy band 0.3 - 8 keV of the field around PSR B2224+65 smoothed with an adaptive 
Gaussian filter. The black cross and white arrow indicate the radio pulsar position and the pulsar's proper motion direction. The 
inset shows the Ha image of the Guitar Nebula as seen by HST WFPC2. 



2. Observations and data analysis 

PSR B2224+65 was observed by the Chandra satellite in 2000 
October 21-22 (Obs ID: 755) with the Advanced CCD Imaging 
Spectrometer (ACIS). The pulsar is located on the back- 
illuminated (BI) ACIS-S3 chip which has a superior quantum 
efficiency among the spectroscopic array. Standard processed 
level-2 data were used in our study. The effective exposure is 
found to be ~ 50 ksec. The frame time of 3.2 sec does not sup- 
port a timing analysis of PSR B2224+65 to search for X-ray 
pulsations. 

The Chandra image of the 4x4 arcmin field around 
PSR B2224+65 is shown in Figure 1. The binning factor of 
the image is 0.5 arcsec. Adaptive smoothing with a Gaussian 
kernel of cr < 2 pixels has been applied to better make vis- 



ible faint diffuse emission. The X-ray point source is clearly 
seen at the radio pulsar position. The inset in Figure 1 shows 
the Ha image observed by the Wide Field Planetary Camera 2 
(WFPC2) on the Hubble Space Telescope (HST). This image 
shows the Guitar Nebula. Romani et al. (1997) reported on a 
4-cr detection of diffuse X-ray emission associated with the 
Ha nebula in ROSAT HRI data. However, we found no evi- 
dence for any diffuse X-ray emission within the correspond- 
ing region of the Chandra image. Rather than an extended X- 
ray feature which has a length of ~ 2 arcmin but which de- 
viates from the pulsar's direction of proper motion by ~ 118° 
is observed. The signal-to-noise ratio of this feature in the en- 
ergy range 2-8 keV is about 6. We further compute the bright- 
ness profile of the feature. From the raw image with bin size 
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Fig. 2. Vignetting corrected brightness profile in the energy 
band 0.3 - 8 keV of the extended feature associated with 
PSR B2224+65. 

Table 1. Proper-motion corrected ephemeris" of 
PSR B2224+65. 



Right Ascension (J2000) 


22' ! 25 m 52.424 s 


Declination (J2000) 


+65° 35' 34.08" 


Pulsar Period, P (s) 


0.682538228276 


Period derivative P (10~ 15 s s -1 ) 


9.55 


Age (10 6 yrs) 


1.13 


Surface dipole magnetic field (10 12 G) 


2.58 


Epoch of Period (MJD) 


49303.00 


DM (pc cnr 3 ) 


36.079 


DM based distance* (kpc) 


1.86 


Spin-down Luminosity (10 33 ) ergs s~ 2 


1.19 



"Manchester et al. (2005). 

b According to Cordes & Lazio (2002). 



of 0.5 arcsec, we estimate counts in eight consecutive boxes 
of 20 x 10 arcsec oriented along the feature. The first box is 
centered at RA=22*25 m 48.609' 5 , Dec=+65°35'45.15" (J2000). 
The vignetting corrected brightness profile in the energy band 
0.3 - 8 keV is shown in Figure 2. The feature was already no- 
ticed by Wong et al. (2003). Inspecting the radio data from the 
NRAO/VLA Sky Survey (NVSS) archive (Condon et al. 1998) 
did not reveal a radio counterpart of it. 

In their brief analysis, Wong et al. (2003) claimed that the 
X-ray emission at the pulsar position is likely to be from a 
shocked nebula. In order to investigate whether there is ex- 
tended nebular emission around PSR B2224+65, we have ex- 
amined the spatial nature of the emission by fitting a con- 
volved 2-D Gaussian to the sub-image around the pulsar po- 
sition with the point spread function (PSF) as the convolu- 
tion kernel. The model PSF at 1 .5 keV is created by using the 
CIAO tool MKPSF. The best-fitted result yields a FWHM of 
0.55 + 0.03 arcsec (l-<x error) which is in agreement with the 



expected width of the Chandra PSF for a point source observed 
on-axis. 

For a spectral analysis of the pulsar emission we extracted 
events from a circle of 3 arcsec radius (encircled energy ~ 
99%) centered on the radio pulsar position. The spectrum from 
the extended X-ray feature was extracted from a box of 1 32 x 1 5 
arcsec, oriented along the direction of this feature. In total, ~ 80 
counts for the pulsar and ~ 370 counts for the feature were se- 
lected in the energy band 0.3 - 8 keV. Even though the photon 
statistics is on the smaller size it is sufficient to provide us with 
interesting insights on the radiation emission nature of the pul- 
sar and the unidentified X-ray feature. 

Response files were computed by using the CIAO tools 
MKRMF and MKARF. According to the photon statistics, the 
spectra were dynamically binned so as to have at least 8 counts 
per bin for the pulsar and at least 30 counts per bin for the ex- 
tended feature. The background spectrum for the spectral anal- 
ysis of the pulsar is extracted from a circle near to the pulsar of 
3 arcsec at RA=22 h 25 m 53.893 s , Dec=+65°35'34.79" (J2000). 
For the analysis of the extended feature, the background spec- 
trum is extracted from a low count region of 18 x 15 arcsec 
centered at RA=22 h 25 m 49.113 s , Dec=+65°36'08.84" (J2000). 
In the energy band 0.3 - 8 keV, the net count rates for the pul- 
sar and the feature are estimated to be (1.53 + 0.18) X 10~ 3 
cts/s and (5.59 + 0.34) x 10~ 3 cts/s respectively. With the aid of 
PIMMS, the pileup fraction of the ACIS-S3 data is estimated 
to be < 0.7% which is negligible. All spectral fits were per- 
formed in the 0.3-8 keV energy range by using XSPEC 11.3.1. 
The degradation of the ACIS quantum efficiency was corrected 
by the XSPEC model ACISABS. The parameters of all fitted 
model spectra are summarized in Table 2. All the quoted errors 
are 1 - <x and were computed for 1 parameter in interest. 

For the pulsar PSR B2224+65 we found that among the 
tested models a single power-law model describes the observed 
spectrum best (y 2 = 0.70 for 6 D.O.F). This model yields a 
column density of Nh < 0.9 x 10 21 crrT 2 , a photon index of 
F = 1.58^ 3 and a normalization at 1 keV of (2.28+°;^) x 10~ 6 
photons keV cm" 2 s _1 . The best-fitted power-law spectrum 
and residuals are shown in Figure 3. 

The unabsorbed fluxes deduced for the best-fitting power- 
law model are fx = 9.3 x 10~ 15 ergs cirT 2 s _1 and fx = 
1.6 x 10~ 14 ergs cirT 2 s _1 in the 0.1-2.4 keV and 0.5-10 keV 
energy bands, respectively. The luminosities in these bands are 
L x =l-l<f 2 kpc x 10 30 ergs s _1 and L x =2.0<i 2 kpc x 10 3() ergs s -1 . 
(iikpc denotes the distance to PSR B2224+65 in units of 1 kpc. 
The conversion efficency, L x (0. 1-2 A keV)/£, is found to be 
0.9<i 2 kpc x 10~ 3 . For a distance close to 1 kpc as suggested by 
Chatterjee & Cordes (2004), this is in good agreement with the 
relation L x (0. 1-2.4 keV)« lO^ found by Becker & Trumper 
(1997) for rotation powered pulsars. 

We have also considered the possibility of a purely ther- 
mal emission scenario. However, fitting a blackbody yields a 
reduced^- 2 of 2.06 for 6 D.O.F. which invalidates this model. 
Although the small number of counts does not support any fit- 
ting with multicomponent models, we were still able to esti- 
mate the upper limit for the polar cap temperature and surface 
temperature by adding a blackbody component to the best- 
fitting power-law model. We computed the confidence con- 
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Table 2. Spectral parameters inferred from fitting the Chandra observed spectrum of PSR B2224+65 and X-ray jet. 



Model 


•) 

Xy 


D.O.F. N H 

10 21 cm" 2 


T/kT 


Normalization h 


PSR B2224+65 


PL 


0.70 


6 


<0.90 


i 5Q+0.43 
JO -0.33 


2.28!^ x 10- 6 


BB 


2.07 


6 


0.00 


0.46 


1.46 x 10- 2 








Extended X-ray 


feature 




PL 


0.79 


9 


< 2.25 


90 +035 


5.66^;^ x 10- 6 


BB 


1.50 


9 


<0.62 


91+ 015 


7.17^-82 x 10 -3 


BREMSS 


0.86 


9 


i 47+I.44 


199 36 +0 ' 64 


2.46:°$ x 10- 5 



" PL = power-law; BB = blackbody; BREMSS = Thermal bremsstrahlung 

* The entry in this column depends on the model in interest. For the power-law model, the unit of the normalization constant is photons keV~' 
cirr 2 s -1 . For the blackbody model, the normalization is (Rkm/D [0 ) 2 , where R km is the source radius in km and D i0 is the distance to the source 
in unit of 10 kpc. For the thermal bremsstrahlung model, the normalization constant is expressed as (3.02 x 10~ 15 /47rD 2 ) where D is 

the source distance in cm and n e and «/ are the electron and ion densities in cm" 3 . 



tours of the blackbody normalization versus the temperature. 
They are shown in Figure 4. Assuming a dipolar magnetic 
field, the conventional size of a polar cap is defined by the last 
open field lines. The radius of the polar cap area is given as 
r pc = R(2nR/cP) l/2 , where R is the neutron star radius, c is 
the speed of light and P is the rotation period of the pulsar. 
For PSR B2224+65, the rotation period of 0.68 s implies a 
polar cap of radius r pc = 175 m. With this estimate we set a 
1 - cr polar cap temperature upper limit of T~ < 1.3 x 10 6 K 
by assuming contribution from one polar cap only. This cor- 
responds to a 1 — cr upper limit on the bolometric luminosity 
of L^j < 6.5 x 10 29 ergs s _1 . Zavlin & Pavlov (2004) have 
estimated a 1 - cr upper limit on the bolometric luminosity of 
Zi*j < 9.1 x 10 29 ergs s _1 by adopting a hydrogen polar cap 
model. However, such models have many free parameters of 
which most are unknown for PSR B2224+65 (e.g. magnetic 
inclination, viewing angle). In view of this, the applicability 
of this model is restricted and a simple blackbody model re- 
garded to provide us with a more conservative upper limit. If 
we assume that the thermal emission is emitted from the whole 
neutron star surface of 10 km radius we compute a surface tem- 
perature upper limit to be T™ < 6.1 x 10 5 K (3 - cr). This limit 
is of the same order as those found for the other old X-ray de- 
tected pulsars (see Becker et al. 2004, 2005, 2006). 

For the linearly extended feature we found that the spec- 
trum is best-fitted by a single power-law model (xl = 0.79 
for 9 D.O.F.). This suggests a non-thermal origin of its emis- 
sion. The model inferred column absorption is found to be 
N H < 2.25 x 10 21 cm 2 . The photon index is T = 0.90+°;^ and 
the normalization at 1 keV is (5.66^ ^) x 10~ 6 photons keV~' 
cm~ 2 s _1 . The best-fitted power-law spectrum and residuals are 
shown in Figure 5. The unabsorbed fluxes of the extended fea- 
ture deduced for the best-fitted model parameters and the en- 
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Fig. 3. Energy spectrum of PSR B2224+65 as observed with 
the Chandra ACIS-S3 detector and fitted to an absorbed power- 
law model (upper panel) and contribution to the x 2 fit statistic 
(lower panel). 



ergy ranges 0.1-2.4 keV and 0.5-10 keV are fx = 2.1 x 10~ 14 
ergs cirr 2 s~' and fx = 1.0 x 10~ 13 ergs cirr 2 s _1 , respec- 
tively. The corresponding luminosities in the 0.1-2.4 keV and 
0.5-10 keV bands are L x = 2.5J 2 kpc X 10 30 ergs s" 1 and 
L x = 1.2af 2 kpc x 10 31 ergs s _1 , respectively. Fitting the spectrum 
of the X-ray feature with the blackbody model does not yield an 
acceptable goodness-of-fit (y 2 = 1.50 for 9 D.O.F.). We have 
also tried to fit the spectrum with a thermal bremsstrahlung 
model. However, the model implies a temperature as high as 
200 keV and hence can be rejected simply because it is not 
physical. 
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Fig. 4. Portion of the l<x, 2<x and 3<x confidence contours 
showing the blackbody normalization versus temperature for a 
blackbody model added on the best-fitted power law model. The 
horizontal line at a normalization of 3.06 corresponds to a polar 
cap radius of 175 m and a pulsar distance of 1 kpc as suggested 
by Chatterjee & Cordes (2004). 
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Fig. 5. Energy spectrum of extended feature detected in the 
field of PSR B2224+65, as observed with Chandra ACIS-S3 
detector and fitted to an absorbed power-law model (upper 
panel) and contribution to the^- 2 fit statistic (lower panel). 

3. Discussion & Conclusion 

According to our analysis, the observed X-rays of 
PSR B2224+65 are of non-thermal origin, suggesting a 
magnetospheric radiation dominated emission scenario. To 
model the detected X-ray spectrum with a blackbody model 
yields no acceptable results. A distinct extended X-ray feature 
of 2 arcmin size is observed to start at the pulsar position. Its 
X-ray spectrum is very hard and no counterpart of it is found in 
Ha nor in NVSS radio data. Comparing the energy spectra of 
PSR B2224+65 and of the extended feature it is obvious that 
the X-ray emission of the later is much harder. We found that 
neither a blackbody nor a thermal bremsstrahlung model yields 
a valid description for its observed spectrum. This contradicts 



to the findings reported by Wong et al. (2003) who fitted the 
emission from the extended feature with a blackbody. That the 
X-ray feature is not aligned with the pulsar's proper motion 
does not a priori invalidate its association with the pulsar. 
Markwardt & Ogelman (1995) reported the ROSAT PSPC ob- 
servation of a highly extended jet-like feature associated with 
the Vela pulsar which also is found to be entirely misaligned 
with the pulsar's proper motion vector. The recent discovery 
of the X-ray trail associated with the millisecond pulsar PSR 
J2124-3358 is another example (Hui & Becker 2006). Such 
asymmetry might be caused by e.g. the inhomogeneities of 
the surrounding ISM or anisotropics in the pulsar wind and/or 
bulk flow of ambient gas (see also Gaensler, Jones & Stappers 
2002). The fact that the feature starts exactly at the position 
of PSR B2224+65 argues against an interpretation as random 
ISM feature. Deeper future observations though are badly 
needed to further study this unique system. 

Same as PSR B2224+65 the X-ray emission of many old 
pulsars are found to be described by power-law spectral mod- 
els (Becker et al. 2004, 2005, 2006) suggesting that their X-ray 
emission is largely dominated by non-thermal emission pro- 
cesses. This is further supported by the fact that the pulse pro- 
files of this class of pulsars are not sinusoidal as it would be 
expected from spin-modulated thermal emission (cf. Becker et 
al. 2004, 2005, 2006). For Bl 133+16 (Kargaltsev et al. 2006) 
and B0943+10 (Zhang et al. 2005), Owing to the limited pho- 
ton statistics, their spectra can be fitted equally well by both 
blackbody and power-law models. While we cannot discrimi- 
nate whether the X-ray emission are thermal or non-thermal in 
these cases, we have not taken these two pulsars in the follow- 
ing discussion. 

According to Cheng & Zhang (1999; hereafter CZ99), the 
non-thermal X-rays from pulsars are synchrotron radiation of 
electron-positron pairs which are created in the strong mag- 
netic field near the neutron star surface by curvature photons. 
These photons are emitted by charged particles on their way 
from the outer-gap to the stellar surface. The fractional size of 
the outer-gap, fo, is defined as the ratio between the mean verti- 
cal separation of the gap boundaries in the plane spanned by the 
rotation axis and the magnetic axis to the light cylinder radius. 
CZ99 estimated the fractional gap size as f = 5.5P 26/21 B^ /? 
and determined the existence of the outer-gap by the criteria 
fo < 1 . We calculated the fo for six recently studied old pul- 
sars and tabulated them in Table 3. Five out of six pulsars have 
fo > 1 which indicates that all but one of these old pulsars 
do not have outer-gaps in the framework of this model. Rather 
then, the CZ99-model predicts that if there is any X-ray emis- 
sion from these old pulsars, it should be thermal. This is ob- 
viously not in agreement with the observations, showing the 
limitation of this model. 

A point, however, which has not been considered in this 
model is that the fractional gap size should depend on the incli- 
nation angle a of the magnetic axis vs. the rotational axis. As 
mentioned by Zhang et al. (2004), the active region of the outer- 
gap should begin at the null charge surface (Q. ■ B = 0) which 
is at the radial distance r,„ from the star. Since the fractional 
size reaches a minimum at r,„, it is more reasonable to deter- 
mine whether the gap exist by the criteria f(r in , a) < 1 rather 
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than by f < I. Especially for larger magnetic inclination an- 
gles the null charge surface is expected to move closer to the 
star, i.e. f(r in ,a) decreases with increasing a. Using Equation 
36 in Zhang et al. (2004) we calculated the variation of /(r,„, a) 
with a for several X-ray detected old pulsars and plotted the re- 
sult in Figure 6. Interestingly, for PSRs B0950+08, B0823+26, 
B0628-28 and B1929+10, the inclination angles deduced from 
fitting the rotating vector model to the radio polarization angle 
swing (Everett & Weisberg 2001; Becker et al. 2005) are all in 
the region /(r,„, a) < 1 (c.f. Figure 6 & Table 3). The emission 
geometry for the PSR B2224+65 is currently unknown which 
prevents us to discuss whether the outer-gap can sustain in this 
pulsar. 

Whether the outer-gap model is of any general meaning for 
the description of the X-ray emission from old pulsars can only 
be constrained by observations. It is necessary to obtain the 
light curves and the spectra from a much larger sample than 
currently available. By now it seems rather speculative to as- 
sume that those X-ray detected old pulsars, by chance, all have 
their inclination angles in that range for which the outer-gap 
model predicts non-thermal emission. 

Acknowledgements. We thank the referee David Helfand for thor- 
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Table 3. The calculated / and f(r in , otrvm) of old pulsars according to CZ99 and Zhang et al. (2004) respectively. 



PSR 


P(s) 


B(10 12 G) 


Q-rvm (degree)" 


fo = 5.5P 2b ' 21 B-; 2 " 


/(^.crvm) 


B2224+65 


0.68 


2.58 




1.99 




B1929+10 


0.23 


0.51 


35.97 


1.31 


0.62 


B0628-28 


1.24 


3.02 
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" The magnetic inclination angles q-rvm are deduced from fitting the rotating vector model (RVM) to the radio polarization angle swing, o-rvm 
of B0628-28 is taken from Becker et al. (2005). Others are adopted from Everett & Weisberg (2001). 
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Fig. 6. Variation of the fractional size of the outer-gap at the null charge surface with magnetic inclination angle for old pulsars. 
The crosses indicate the calculated f(r in , ffRvivi) for the pulsars with measured magnetic inclination angle (cf. Table 3). 



1 I 1 I 1 I 1 I 1 I 1 


1 1 1 1 






- B2224+65 






B1 929+10 






B0628-28 






B0950+08 






B0823+26 






J2043+2740 










^^^^ 






1 . 1 . 1 . 1 . 1 


i,i, 





